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Abstract 

Weak gravitational lensing of distant galaxies by foreground structures has proven 
to be a powerful tool to study the mass distribution in the universe. Nowadays, 
attention has shifted from clusters of galaxies to the statistical properties of the large 
scale structures and the halos of (field) galaxies. These applications have become 
feasible with the advent of panoramic cameras on 4m class telescopes. 

In this review we will give an overview of recent advances in this fast evolving 
field of astronomy. We start with a discussion of the recent measurements of weak 
lensing by large scale structure ("cosmic shear"), which can be used to constrain 
cosmological parameters. We also show how weak lensing can be used to measure 
the relation between galaxies and dark matter (galaxy biasing) directly. 

Another area that benefitted greatly from the current data is the weak lensing 
by galaxies (galaxy-galaxy lensing). Weak lensing provides a unique probe of the 
gravitational potential on large scales. Hence, in the context of dark matter, it can 
provide constrains on the extent and shapes of dark matter halos. Furthermore, it 
can test alternative theories of gravity (without dark matter). The best constraint 
comes from the first detection of the anisotropy of the lensing signal around lens 
galaxies, which suggest that the dark matter halos are flattened. An isotropic signal, 
as predicted by MOND, is excluded with 99% confidence. 
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1 Introduction 



Weak gravitational lensing is a valuable tool to study the (dark) matter dis- 
tribution of a range of objects in the universe. The differential deflection of 
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light rays by intervening structures allows us to study the projected mass dis- 
tribution of the deflectors, without having to rely on assumptions about the 
state or nature of the deflecting matter. 



In 1984 Tyson et al. tried to detect the weak lensing signal induced by an 
ensemble of galaxies using photographic plates. However, the first succesful 
detection was made using CCD images of the massive cluster of galaxies A1689 
( [Tyson et al., 1990|) . Although much theoretical work concentrated on weak 
lensing by large scale structure (e.g., |lVliraida-Escude, 199 1|; [Blandtord et al 



1991| ; [Kaiser, 1992| ), this area of astronomy blossomed with study of the mass 
distribution in clusters of galaxies (e.g., [Bonnet et al., 1994[ ; [Fahlman et al. 



1994[ ; p^uppino fc Kaiser, 1997[ ; pioekstra et al., 1998[ , pOO0| , for an extensive 
review see Mellier 1999) and galaxy groups ([Hoekstra et al., 2001a[) . These 
results were an important first step in demonstrating the feasibility of the 
technique, but nowadays more and more studies concentrate on the "field". 

Adequate schemes to correct for a variety of observational distortions (such as 
PSF anisotropy, seeing, camera distortion) were developed during the 1990s 
(e.g.. Kaiser et al., 1995 ; Luppino fc Kaiser, 1997[ ; Hoekstra et al., 1998 ; Kui- 
jken, 1999[ ; [Kaiser, 2000| ; Pernstein fc Jarvis, 2002[) . Together with the de- 
velopment of wide field CCD cameras, these advances allowed for the first 



detections of weak lensing by large scale structure in the spring of 2000 ([Ba- 
con et al., 2000| ; [Kaiser et al., 2000[ ; [Van Waerbeke et al, 2000[ ; [Wittman et 



"20001) . 



These first detections were based on small data sets, and consequently the 
accuracy of the measurements was low. However, the amount of data is in- 



creasing rapidly (e.g.. Bacon et al., 2002; 


Maoli et al., 2001; Hoekstra et al.. 


2002a 


,[b[; [Refregier et al., 2002 


; Van Waerbeke et al, 2001a 


. .'(KI J 


). Current data 



sets cover a few to tens of degrees on the sky. Planned surveys, such as the 
CFHT Legacy Survey, which will obtain deep imaging data for 170 deg^ in 
five filters, will be another major leap forward. The constraints on cosmologi- 
cal parameters are tightening rapidly, and the consistency of the most recent 
results ( [Bacon et al., 2002[ ; [Hoekstra et al., 2002 b[ ; [Refregier et al., 2002| ; [Van 
Waerbeke et al, 2002|) suggest that cosmic shear can play an important role 
in this era of "precision cosmology" . 



Another important application of weak lensing is the study of the distribu- 
tion of the dark matter with respect to the galaxies. For instance, one can 
study the statistical properties of the galaxy and dark matter distributions: 
the biasing relations ( [Hoekstra et al., 2001b[ , [2002c[ ) . Such measurements can 
provide important constraints on models of galaxy formation. Weak lensing 
can also be used to study the properties of the dark matter halos surrounding 
galaxies. Rotation curves of spiral galaxies have provided important evidence 
for the existence of dark matter halos (e.g., [van Albada fc Sancisi, 195B[). Also, 
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strong lensing studies of multiple imaged systems require massive halos to ex- 
plain the oberved image separations. However, both methods provide mainly 
constraints on the halo properties at relatively small radii. 



The weak lensing signal can be measured out to large projected distances, 
and in principle it can be a powerful probe of the potential at large radii. 



constraining the amount of dark matter (e.g., Brainerd et al., 1996; Grithths 


et al., 1996; Dell'Antonio & Tyson, 199q; |Hudson et al., 1998 


Fischer et al., 


200q; tloekstra, 2000|; |Wilson et al., 2001; McKay et al., 2001 


, the extent of 


the dark matter halos (e.g.. 


Brainerd ct al., 1996; Hudson et al., 1998|; [Fischei 


et al., 200q; [Hoekstra, 200^; 


Hoekstra et al., 2002d,d|) or the a^ 


Aerage shape of 


the halos (|Hoekstra et al., 2002d). 



So far, we have assumed that we can apply the theory of General Relativity 
to describe gravity. Consequently, the observational data require the existence 
of dark matter. However, at the large scales probed in this paper, alternative 
theories have been proposed to explain the shapes of galaxy rotation curves 
without the use of dark matter. In particular Modified Newtonian Dynamics 
(MOND; IMilgrom (1983| )) succesfuUy reprod uces rotation curves using only 
visible matter. Ideally, one would like to test such alternative theories far away 
from the visible matter. Potentially, weak lensing can place tight constraints 
on such theories ( [Mortlock fc I'urner, 2001a| JBD. 

In this review, we will address various applications of, and recent developments 
in weak gravitational lensing. In Section 2 we concentrate on the cosmic shear 
measurements, and Section 3 deals with galaxy biasing. The constraints on the 
halos of galaxies from weak lensing are described in Section 4. In Section 5 we 
discuss how the lensing data can be used to place limits on alternative theories 
of gravity. 

Throughout the paper we will make use of results based on i^p-band imaging 
data from the Red-Sequence Cluster Survey (RCS; IGladders fc Yee, 2000| ; 
Yee fc Gladders, 2001| ). We refer the reader to [Hoekstra et al. (2002a|) for a 
thorough discussion of the weak lensing analysis of these data. 



2 Lensing by large scale structure 

The measurement of the coherent distortions of the images of faint galaxies 
caused by weak lensing by intervening large scale structures provides a di- 
rect way to study the statistical properties of the growth of structure in the 
universe. Compared to many other methods, weak lensing probes the matter 
directly, regardless of the light distribution. In addition it provides measure- 
ments on scales from the quasi-linear to the non-linear regime where compar- 
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isons between observations and predictions are still limited. 



We provide a basic description of the theory of weak lensing by large scale 
structure. Detailed discussions on this subject can be found elsewhere (e.g., 
Schneider et al., 199^ ; Bartelmann fc Schneider, 2001 ). Current measurements 
of cosmic shear have concentrated on the two-point statistics, which can be 
related to the convergence power spectrum, which is defined as 



9Hq 




Jk{w) 



w 



where w is the radial (comoving) coordinate, wh corresponds to the horizon, 
a{w) the cosmic scale factor, and fxiw) the comoving angular diameter dis- 
tance. As shown by [Jain fc Seljak (1997] ) and [Schneider et al. (1998] ) it is 
necessary to use the non-linear power spectrum in equation (1). This power 
spectrum is derived from the linear power spectrum following the prescrip- 
tions from [Peacock fc Dodds (1996|) . We note that this approach might not be 
accurate enough for future measurements (e.g., [Van Waerbeke et al, 200 lb| , 
200^ ). 



W{w) is the source-averaged ratio of angular diameter distances Dis/Ds for a 
redshift distribution of sources Pb{w): 



W{w)= / dwpbiw) — . . 



(2) 



Hence, it is important to know the redshift distribution of the sources, in order 
to relate the observed lensing signal to Pk(0- 

Until recently, the top-hat variance was the most widely used two-point statis- 
tic. It is related to the convergence power spectrum through 
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(3) 



where 6 is the radius of the aperture used to compute the variance, and Ji is 
the first Bessel function of the first kind. 



Another useful statistic is the aperture mass, which is defined as (e.g., [Kaiser 
et al., 1994[ ; [Schneider et al., 1998| ) 



M,^{e) = / d'<pU{(t>)^^{<P). 



(4) 
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Provided f/(0) is a compensated filter, i.e., / d(t)(l)U {(p) = 0, with = for 
4> > 0, the aperture mass can be expressed in term of the observable tangential 
shear 7t using a different filter function Q{(j)) (which is a function of U{(f))). 
Using the choice of filter functions from Schneider et al. (1998| ), the variance 
of the aperture mass (M^p) is related to the power spectrum through 



CM 

(M,y = 2n J dl IP^{1) 



12 

Ti{ie) 



(5) 



where J4 is the fourth-order Bessel function of the first kind. The measurement 
of (M^p) is of particular interest, because the measurements at different scales 
are only weakly correlated. 



The first detections of cosmic shear ( |Bacon et al., 200q ; [Kaiser et al., 200^ ; 
IMaoh et al., 200 1| ; [Van Waerbeke et al, 2000| ; [Wittman et al., 200q) actually 



computed the excess variance in apertures directly from the data. Some more 
recent studies have used the same technique ( [Bacon et al ., 2002|; [Hoekstra et 



al., 2002a| ; [Refregier et al, 2002| ; [Van Waerbeke et al, 2001aj) 



However, recent studies show that an optimal use of the data is to measure the 
shear correlation functions from the data ([Pen et al., 20 02 



Van Waerbeke et 



al, 2002| ; [Hoekstra et al., 2002 b[ ). These correlation functions can be related to 
the various two-point statistics ([Crittenden et al., 2002| ; Pen et al., 2002[) . In 
addition, this approach allows one to split the signal into two components: an 
"E"-mode, which is curl-free, and a "B"-mode, which is sensitive to the curl of 
the shear field. Gravitational lensing arises from a gravitational potential, and 
hence it is expected to produce a curl-free shear field. Thus, the "B"-mode can 
be used to quantify the level of systematics involved in the measurement. The 
decomposition is naturally carried out by using the aperture mass statistic 

Map. 

Several sources of B-mode have been identified. For instance, simple models 
describing the intrinsic alignments of galaxies predict a small B-mode (e.g., 
Crittenden et al., 2002| ), although the amplitude is still uncertain. Hence, 



any measured B-mode is dominated by residual systematics in the data (e.g., 
imperfect correction of the PSF anisotropy) or intrinsic alignments. The effect 
of intrinsic aligments can be minimized by selecting galaxies with a broad 
redshift distribution. In future surveys, with photometric redshift information 
for the galaxies, the contribution from intrinsic aligments can be removed 
completely by correlating the shapes of galaxies with different redshifts. 

Figure 1 shows the results of such a decomposition into "E" and "B" -modes. 
The analysis, which is described in detail in [Hoekstra et al. (2002b[) , uses 53 
deg^ of RCS imaging data. For all three samples we find that the "B"-mode 
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Fig. 1. The upper panels show the measured variance of the aperture mass (Mfp) 
as a function of aperture size 9 for different samples of source galaxies. This signal 
corresponds to the E-mode. The lower panels show the variance {MJ_) when the 
phase of the shear is increased by 7r/2, and corresponds to the B-mode. The error 
bars indicate the la statistical uncertainty in the measurements, and have been 
derived from the field-to- field variation of the 13 RCS patches (thus the error bars 
include cosmic variance). Note that the points are slightly correlated. We detect a 
significant B-mode on scales 5—10 arcminutes. On scales larger than 10 arcminutes 
the B-mode vanishes. The sample of bright galaxies (20 < Rc < 22) should not 
be affected significantly by systematics because their sizes are large compared to 
the PSF. Therefore the significant B-mode at scales of a few arcminutes is likely 
to be caused by intrinsic alingments. To minimize the effect of intrinsic alignments 
on our cosmological parameter estimation, we will use the sample of galaxies with 
22 < Rc < 24 to this end. 

(lower panels) vanishes on scales larger than ~ 10 arcminutes, suggesting 
that neither observational distortions or intrinsic alignments of sources have 
corrupted our measurements. 



On smaller scales we detect a significant "B"-mode. Interestingly, we detect a 
significant signal for the bright galaxies (20 < Rc < 22). These galaxies have 

sizes that are large compared to the PSF, and therefore they are less affected 
by residual systematics. Intrinsic alignments are expected to be particularly 
important for these bright galaxies, and we therefore conclude that (at least 
part of) the observed "B"-mode is likely to be caused by intrinsic alignments. 
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Fig. 2. Left panel: The measured variance of the aperture mass (M^p) as a func- 
tion of aperture size for source galaxies with 22 < Rc < 24. The error bars 
have been increased to account for the unknown correction for the "B"-mode ob- 
served in Figure 1. Right panel: Contraints on Qrn and ag using priors Gaussian 
priors on F, J7tot = + ^^A; and the source redshift distribution Zg (which gives 
{z) = 0.59 lb 0.02). The hkehhood contours have been derived by comparing the 
measurements to CDM models with n = 1. For F we used the constraints from 
the 2dF survey F = 0.21 ± 0.03 ( [Peacock et al., 200 Ij [Efstathiou et al., 20011) , and 
for r^tot we used the Boomerang constraints Otot = 1-02 it 0.06 ( [Netterfield et a~ 
2001 ). The contours indicate the 68.3%, 95.4%, and 99.7% confidence limits on two 



paramaters jointly. 



The next step is to relate the observed signal to estimates of cosmological pa- 
rameters. Our current understanding of the redshift distribution of the source 
galaxies is sufficient to obtain accurate constraints. The most accurate re- 
sults are derived when one uses external priors on the parameters. Useful 
constraints come from redshift surveys and the temperature fluctuations in 
the CMB. Current weak lensing measurements provide joint constraints on 
flrn and cTg, and Figure 2 shows the results from Poekstra et al. (2002bD . How- 
ever, measurements of higher order statistics, such as the skewness, can be 
used to break the degeneracies between flm, and as (e.g., [Van Waerbeke et| 
al., 19991) . Recently the first detection of a measure of skewness was reported 
([Bernardeau et al., 20024 |b|), and the prospects for measurements of higher 



order statistics are good. 

In Table 1 we list the estimates of as (for Qm = 0.3 and F = 0.21) as derived 
by the most recent weak lensing surveys ( Paeon et al., 2002| ; Poekstra et 
al., 2002b| ; [Refregier et al., 2002| ; |Van Waerbeke et al, 2002]) . Despite the 
fact that these results have been obtained using a wide variety of telescopes 
(both ground based and HST), filters, and integration times, the agreement is 
remarkable. 



We note that the value of as determined from the RCS data is essentially 
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Cosmic shear survey 



0-8 



Hoekstra et al. (2002b|) (RCS) 0.86+° °^ 



Bacon et al. (2002|) 
Refregier et al. (2002|) 
Van Waerbeke et al (200^) 



0.05 

n 07+0.10 

-0.09 

0.94 ±0.14 
0.98 ±0.06 



Table 1 

Values of ug and 68% confidence intervals as derived from 4 independent cosmic 
shear measurements (adopting 0.^ = 0.3, = 0.7, and T = 0.21. 

constrained by the measurements at scales larger than 10 arcminutes, where 
the "B"-mode is negligible. The result from Van Waerbeke et al. (2002) might 
be biased high, because of their large scale "B"-mode. Both [Bacon et al 



(2002|) and [Refregier et al. (2002|) do not separate their signal into "E" and 
"B '-modes, and therefore it might also include some residual systematics. 



A widely used method to determine the normalization of the power spectrum 
uses the number density of rich clusters of galaxies (e.g., [Borgani et al., 2001 



Carlberg et al., 1997|; |Eke et al., 1996|; Fan & Bahcall, 1998; Pen, 1998a; 


Pierpaoli et al., 2001; 


Reiprich & Bohringer, 2001; 


Seljak, 2001; jViana et al., 



200 1| ). Such systems are rare, and as a result a very sensitive probe of as, 



provided one can determine their mass. The derived values for as from this 
technique range from values as low as 0.61 ±0.05 ( [Viana et al., 200"l| ) to values 
around unity (e.g., [Fan fc Bahcall, 1998| ; pen, 1998a| ; [Pierpaoli et al., 2001[) . 
The statistical error bars on most of these measurements are small, and hence 
the large spread in values suggests an underlying large systematic uncertainty. 



3 Galaxy biasing 



The growth of structures in the universe via gravitational instability is an 
important ingredient in our understanding of galaxy formation. However, the 
connection to observations is not straightforward, as we need to understand 
the relation between the dark matter distribution and the galaxies themselves. 
Galaxy formation is a complex process, and it is not guaranteed a priori that 
this relation, referred to as galaxy biasing, is a simple one. The bias might be 
non-linear, scale dependent or stochastic. In the simplest case, linear, deter- 
ministic biasing, the relation between the dark matter and the galaxies can 
be characterized by a single number b (e.g., [Kaiser, 1987[ ). 



Most observational constraints of biasing come from dynamical studies (see 
Strauss fc Willick, 1995[) which probe relatively large scales (lO/i^g^ Mpc or 



more). Recent estimates on these scales suggest values of 6 ~ 1 for galaxies 
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e.g., Peacock et al., 2001| ; |Verde et al., 20011) . On smaller scales some con- 



straints come from measurements of the galaxy two-point correlation function, 
which is compared to the (dark) matter correlation function computed from 
numerical simulations. These studies indicate that the bias parameter 6 ~ 0.7 
on scales less than ~ 2h'^Q Mpc (Vim = 0.3, = 0.7). On larger scales b 
increases to a value close to unity ([Jenkins et al., 1998|) Although this proce- 
dure provides useful information about the bias parameter 6 as a function of 
scale, it does rely on the assumptions made for the numerical simulations. In 
addition, it cannot be used to examine how tight the correlation between the 
matter and light distribution is. To do so, we need to measure the galaxy-mass 
cross-correlation coefficient r, which is a measure of the amount of stochastic 
and non-linear biasing (e.g., [Pen, 1998b| ; pekel fc Lahav, 1999| ; ^omerville et 



al., 2001[) . 



Redshift surveys can be used to determine the relative values of b and r for 
different galaxy types ( [legmark fc Bromley, 19"99| ; [Blanton, 2000| ), but weak 



gravitational lensing provides the only direct way to measure the galaxy-mass 
cross-correlation function (e.g., [l''ischer et al., 2000| ; [Wilson et al., 200T[ ; [McKay 



et al., 200 1[ ; [Hoekstra et al., 200 lb[) . [Fischer et al. (200q ) used the SDSS 



commissioning data to measure the galaxy-mass correlation function and their 
results suggested an average value of 6/r ~ 1 on submegaparsec scales. This 
approach has been explored by [Guzik fc Seljak (200 1[ ) who used semi-analytic 
models of galaxy formation combined with N-body simulations. Their results 
suggest that the cross-correlation coefficient is close to unity. 

To study the galaxy biasing we use a combination of the galaxy and mass 
auto-correlation functions, as well as the cross-correlation function. A com- 
mon definition of the "bias" parameter is the ratio of the variances of the 
galaxy and dark matter densities, which is the definition we will use here. In 
the case of deterministic, linear biasing the galaxy density contrast 6g is simply 
related to the mass density contrast 6 as 6g = b5 ( [Kaiser, 1987[ ), and the ratio 



of the variances is the only relevant parameter. However, the bias relation is 
likely to be more complicated as it depends on the process of galaxy formation, 
and might be stochastic, non-linear or both. We allow for non-linear stochastic 
biasing by including the galaxy-mass cross-correlation coefficient r (e.g., [Fen 



1998b[ ; Pekel fc Lahav, 1999[ ; ^omerville et al., 200 1[ ). The cross-correlation 



coefficient r mixes non-linear and stochastic effects ( Pekel fc Lahav, 1999| ; 



Somerville et al., 2001[ ) which we currently cannot disentangle with weak lens- 



mg. 



The bias parameters are related to the observed correlation functions through 
([Hoekstra et al., 2002c[) 



/l(gap,^^^,^^A)x (6) 
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and 



where /i and /2 depend on the assumed cosmological model and the redshift 
distributions of the lenses and the sources (see [Van Waerbeke, 1998| ; [Hoekstra 



et al., 2002c|). The values of /i and /2, however, depend minimally on the 



assumed power spectrum and the angular scale. 

There is no reason for 6 or r to be constant with scale, but as long as h and 
r vary slowly with scale, we can still infer the bias parameters from Eqs. (3) 
and (4). The Map statistic is sensitive to a fairly small range in k in Fourier 
space for a given aperture size. If the galaxy power spectrum Pg^{k) can be 
related to the matter power spectrum Pmm(^) through P^{k) = b'^{k)P^m{k), 
we essentially measure the average value of b{k) in the A;-range probed by the 
Map statistic ( [Van Waerbeke, 1998| ; [Hoekstra et al., 2002c[) . In other words. 



Map is a pass-band filter, which also explains why different scales are only 
slightly correlated ( [Schneider et al., 1998[ ). 



Hoekstra et al. (2002c) studied the bias properties of galaxies with 19.5 < 
Rc < 21. They combined the weak lensing measurements from the RCS and 
the VIRMOS-DESCART survey (e.g., [Van Waerbeke et al, 2002[) to obtain the 
first direct measurements of the bias parameter b and the galaxy-mass cross- 
correlation coefficient r as a function of scale (with scales ranging from 0.1 to 
4 Mpc). The results of this analysis are presented in Figure 3. [Hoekstra et 



al. (2002c[ ) find that both b and r vary with scale for this sample of galaxies. 



The low value of r on scales ~ 0.5 — lh~^ Mpc suggests significant stochastic 
biasing and/or non-linear biasing. Although both b and r vary with scale, the 
ratio b/r turns out to be almost constant with scale. Hoekstra et al. (2002c[) 
find that b/r = 1.090 ± 0.035 out to Yh'^ Mpc. 

Models of galaxy formation predict the relation between the galaxies and the 
underlying dark matter distribution. The two commonly used approaches are 
hydrodynamic simulations (e.g., Blanton et al., 20DI1| ; Yoshikawa et al., 2001) or 



e.g. , [Kauft'mann et al., 1999a[ ,[E[; [Somerville et al., 200 1[ ; [Guzik fc Seljak, 2001 



A comparison of the models with the lensing results can provide accurate 
constraints. Unfortunately, current models of galaxy formation are limited to 
relatively bright (i.e., massive) galaxies, because of the mass resolution of the 
numerical simulations. Most of the galaxies in the sample studied by [Hoekstra 



et al. (2002c[) are too faint to be included in the simulations. This severely 



limits a direct comparison. However, with available multi-color data the weak 
lensing results can be matched better to the models, and vice versa, future 
models should include lower mass galaxies. 
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Fig. 3. (a) The measured value of the galaxy-mass cross correlation coefficient r 
as a function of scale for the ACDM cosmology, (b) The bias parameter 6 as a 
function of scale. The upper axis indicates the effective physical scale probed by 
the compensated filter at the median redshift of the lenses {z = 0.35). The error 
bars correspond to the 68% confidence intervals. Note that the measurements at 
different scales are slightly correlated. 
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4 Constraints on the halos of galaxies 



Rotation curves of spiral galaxies (e.g., |van Albada fc Sancisi, 19861 ) and strong 



lensing analyses of multiple imaged sources have provided strong evidence for 
the existence of massive dark matter halos around galaxies. However, the 
lack of visible tracers hampers our knowledge of the gravitational potential 
at large projected distances. Only satellite galaxies (e.g., |Zaritsky fc White] 
1994| ) provide a way to probe the outskirts of isolated galaxy halos. 



However, the weak gravitational lensing signal induced by such galaxies can 
be measured out to large projected distances, thus providing a powerful probe 
of the potential at large radii, which can be used to study the extent of the 
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dark halos (e.g., Prainerd et al., 1996| Pudson et al., 1998| ; [Fischer et al., 2000| ; 



Hoekstra, 2000| ; [Hoekstra et al., 2002d| ,|dl). 



The lensing signal induced by an individual galaxy is too low to be detected, 
and one has to study the ensemble averaged signal around a large number of 
lenses. Early measurements of the galaxy-galaxy lensing signal were limited by 
the small number of lenses/sources (e.g., Prainerd et al., 1996| ; Pudson et al.J 



199(j| ), because of the lack of panoramic cameras. Recently the accuracy with 
which the lensing signal can be measured has improved significantly. [Fischer et 



al. (200(][ ) measured a very significant signal, and [Wilson et al. (20(JI| ) studied 
the signal around early type galaxies as a function of redshift. [McKay et al.[ 
(2001[) used the redshift information from the SDSS to study the galaxy-galaxy 
lensing signal function of galaxy properties. 

Here we present some preliminary results based on 45 deg^ of RCS Rc band 
data. The detailed analysis is described in [Hoekstra et al. (2002d[ ), and a 



description of the method can be found in Poekstra (20001) . We use galaxies 



with 19.5 < Rc < 21 as lenses, and galaxies with 21.5 < Rc < 24 as sources 
which are used to measure the lensing signal. The simplest approach to galaxy- 
galaxy lensing is to measure the ensemble averaged tangential distortion as 
a function of radius around the sample of lenses. The result for our data is 
presented in Figure 4a. The signal when the phase of the shear is increased by 
7r/2 is presented in Figure 4b. This signal is consistent with zero, indicating 
that observational systematics have been successfully corrected for. 

We fit a SIS model to the tangential distortion which yields {te) = 0'.'129 ± 
O'.'Oll. Given the redshifts of the lenses and the sources this corresponds to 
a velocity dispersion of (cr^)^/^ = 123 ± 5 km/s. The corresponding circular 
velocity can be obtained using Vc = V2a. 

The tangential shear profile is a convolution of the actual halo profile around 
galaxies and the clustering properties of the lenses. Pischer et al. (2000| ) used 



this tangential shear profile, taking into account the clustering of the lenses, 
to obtain constraints on the extent of the dark matter halos surrounding the 
lenses. They conclude that the halos are large, but they were unable to obtain 
tight constraints. In Section 4.1 we show the results of a maximum likelihood 
analysis of the data, and demonstrate that this approach gives much better 
constraints on the sizes of halos. 

An application of galaxy-galaxy lensing that is still in its early stage is the 
measurement of the shapes of galaxy halos. If the halos are flattened, the 
lensing signal will be slightly anisotropic. Provided the halos and the galaxies 
are aligned, we can measure this angular dependence of the tangential shear. 
In Section 4.2, we present the first results of such an analysis. 
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Fig. 4. (a) The ensemble averaged tangential distortion as a function of radius 
around the sample of lens galaxies. The measurements have been corrected for the 
presence of source galaxies associated with the lens galaxies, (b) The signal when 
the phase of the distortion is increased by 7r/2: no signal should be present if the 
signal in (a) is due to lensing. The vertical error bars indicate the la errors 

4-1 Sizes of galaxy halos 



For an isolated lens, the induced lensing signal is purely tangential. This is no 
longer the case for an ensemble of lenses (which also cluster) : the other lenses 
cause small perturbations, and the shear is no longer tangential. To study the 
sizes of halos around galaxies [l^'ischer et al. (20001) used only the tangential 



component. However, Hoekstra (2000| ) was able to obtain better constraints 



using a maximum likelihood analysis on their imaging data of CN0C2 fields. 
The maximum likelihood analysis uses a parameterized model for the mass 
distribution of individual galaxies is compared to the data. This has the ad- 
vantage that one uses the information contained in both components of the 
distortion. It implicitely assumes that all clustered dark matter is associated 
with galaxies. Furthermore, the contributions from faint neighboring galaxies 
are not included in the model, and as such, the results reflect the average 
density profile around the lenses. 



A useful model to describe a truncated halo is ([Schneider fc Rix, 1997| ; [Hoek- 
stra, 200q ) 



2Gr \ + 



where s is a measure of the truncation radius (i.e., the scale where the pro- 
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Fig. 5. Likelihood contours for the velocity dispersion a^, and truncation parameter 
of a galaxy with a fiducial luminosity = IO^'^/i^^Lb©- To derive this result 
we assumed that the velocity dispersion a scales oc L^^, and that the truncation 

parameter s scales oc L^"^ (or oc cr^). The contours indicate the 68.3%, 95.4%, and 
the 99.7% confidence on two parameters jointly. 

file steepens). The total mass of this model is finite, and half of the mass is 
contained with r 



As before we assume that a oc L^^^. We also have to assume a scaling relation 
for the truncation parameter s. Currently, no such observational constraints 
exist, although galaxy-galaxy lensing analyses with redshifts for the lenses can 
be used to this end. As shown by Poekstra (2000| ) different scaling relations 
lead to somewhat different values for s. Also one can compare the predictions 
from N-body simulations ([Navarro et al., 1997| ) to the data. The aim here is 



not to derive the definitive value of s, but to demonstrate that weak lensing 
can constrain the "extent" of the dark matter halos. 

Currently, we do not have (photometric) redshift information for the lens 
galaxies in the RCS fields. To derive constraints on the dark matter halos we 
follow bioekstra (20001 ), and make mock catalogs of lenses where the redshifts 
are drawn from the CN0C2 redshift survey ([Vee et al., 2000|) on the basis of 
their Rc-hand magnitude. This provides a very crude redshift estimate. We 
found that the derived results did not change much between realisations. 

In the analysis we assumed that the velocity dispersion a scales oc L^^, and 
that the truncation parameter s scales oc L^^ (or oc a^). Figure 5 shows the 
results of maximum likelihood analysis for a galaxy with a fiducial luminosity 
of = 1O^°/i~^Lb0. For such a galaxy we find a* = 118t4, f^iir agreement 
with other estimates. In addition, the average extent of the dark matter halo 
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is well constrained, and we obtain s* = 265125^ ^ kpc, and a 99.7% confidence 
upper limit of < 390h~^ kpc. 



4-2 Shapes of galaxy halos 



Numerical simulations of cold dark matter result in triaxial halos, with typical 
ellipticity of ~ 0.3 (e.g., Pubinski fc Carlberg, 19911) . In the context of colli- 



sionless cold dark matter, the theoretical evidence for flattened halos is quite 
strong. However, the observational evidence is still limited, ^ackett (1999| ) 
gives an overview of current constraints which probe the vertical potential on 
scale < 15 kpc. These results suggest an average value of c/a = 0.5 ± 0.2 
(where c/a is the ratio of the shortest to longest principle axis of the halo). 

Current constraints are mainly limited by the lack of visible tracers that can 
probe the gravitational potential around galaxies. Hence, galaxy-galaxy lens- 
ing is potentially the most powerful way to derive constraints on the average 
shape of dark matter halos. As such, it provides an important test of the colli- 
sionless cold dark matter paradigm. [Brainerd fc Wright (2000| ) and |JNatarajan 



fc Refregier (20001) examined the data requirements needed for a detection of 
the halo flattening. Although their estimates appear to be too optimistic, it 
is clear that large amounts of data are required. The measurement of the av- 
erage shape of the dark matter halos is much more difficult than the previous 
measurements: the galaxy-galaxy lensing signal is small, and now one needs to 
measure an even smaller azimuthal variation. Furthermore one has to assume 
that the galaxy and its halo are aligned, which appears to be a reasonable 
assumption. 

Our approach is different from that of [Brainerd fc Wright (2000|) and [Natara- 



jan fc Refregier (2000 ), who propose to study the azimuthal variation in the 
tangential shear around the lenses. Instead we use a maximum likelihood anal- 
ysis to derive constraints on the shapes of dark matter halos. We note that 
this is a first attempt to do so, and more detailed work is required. We simply 
assume that the (projected) ellipticity of the dark matter halo is related to the 
shape of the galaxy as Chaio = f^iens, and use the best fit parameters for the 
velocity dispersion and truncation. The ellipticities of the lenses have been cor- 
rected for the PSF, similar to the sources. The lensing signal is dominated by 
the contribution from early type galaxies, and therefore our parameterisation 
of the halo shape is reasonable. 

The resulting Ax^ as a function of / is presented in Figure 6. The main 
conclusion from this analysis is that spherical halos are excluded with 99% 
confidence. Although the inferred ellipticity of dark matter halos depends on 
the adopted dependence on the light distribution, we note that the exclusion 
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Fig. 6. Ax^ as a function of /. We have assumed that the eUipticity of the halos is 
related to the observed elhpticity of the lens as e^aio = /eicns- We have indicated the 
68.3% and 95.4% confidence intervals. Round halos (/ = 0) are excluded with 99% 
confidence. The average elhpticity of the galaxies is (eiens) = 0.261 and hence the 
average halo elhpticity is (ehaio) = 0.2llQ Qg (68.3% confidence), which corresponds 
to an average projected axis ratio of c/a = 0.65^002 (68.3% confidence). 

of a spherical halo is robust. Compared to the measurement of the tangen- 
tial shear as a function of radius, the result presented here is more sensitive 
to residual systematics. A detailed analysis, however, demonstrates that our 
result is robust, and not caused by systematics. 

The average eUipticity of the lens galaxies is (eiens) = 0.261. The best fit value 
of / = 0.8 then implies a average projected halo elhpticity of (chaio) = 0.21i'^q;q8 
(68.3% confidence), which corresponds to an projected axis ratio of c/a = 
0.65lo;o2 (68.3% confidence). Although the weak lensing yields a projected 
axis ratio, the result is in fair agreement with the numerical simulations. We 
note, however, that it is not straightforward to interpret this result: the lensing 
signal comes from a range of galaxy types, although typically the lensing signal 
is dominated by the early type galaxies, for which our halo model might be 
appropriate. 



5 Limits on alternative theories of gravity 

In the context of General Relativity, we have to invoke large amounts of dark 
matter to explain a wide variety of observational results. However, at the large 
scales probed in this paper, no definite tests of general relativity have been 
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made, and alternative theories of gravity have been proposed to explain the 
shapes of galaxy rotation curves without the use of dark matter. In particular 
Modified Newtonian Dynamics (MOND; [Milgrom (1983| )) succesfully repro- 
duces rotation curves using only visible matter. 

Dynamical studies can provide useful constraints on such theories, but only 
in the regions where there is visible matter. The best regime to test such 
alternative theories is far away from the visible matter. This requires a way to 
probe the gravitational potential on large physical scales. Weak gravitational 
lensing is currently the best approach to this end. As shown above, it can 
accurately probe the gravitational potential on scales where other methods 
fail. In addition, in the absence of dark matter halos, galaxy-galaxy lensing 
is simple, because the galaxies can be treated as point masses. Unfortunately 
most alternative theories of gravity do not provide a useful description of 
gravitational lensing. For instance, MOND lacks a relativistic theory that can 
describe hght deflection. 

To test alternative theories of gravity we can use two different approaches. 
A measurement of the radial dependence of the lensing signal gives the best 
accuracy, but it requires knowledge of the deflection law. The most direct test 
is the detection of the the azimuthal variation of the lensing signal around 
the lenses: alternative theories predict an isotropic lensing signal, and the 
measurement does not require knowledge of the deflection law! 

In alternative theories of gravity, any anisotropy in the lensing signal caused 
by the intrinsic shapes of the galaxies (i.e., the light and gas distribition) 
decreases oc r~^, and hence is negligible on the scales probed here. Therefore 
the observed anisotropy in the galaxy-galaxy lensing signal (see Section 4.2) 
poses a serious problem for alternative theories of gravity. More detailed work 
is required to confirm this result, but new data sets will provide much better 
constraints in the very near future, and allow us to study the anisotropy a 
function of projected distance from the galaxy. 

The detection of the anisotropy of the lensing signal around galaxies provides 
the best evidence for the existence of dark matter halos around galaxies. For 
completeness we also discuss constraints on the radial dependence of gravity, 
following the approach suggested by [Mortlock fc Turner (2001b| ) . 

5. 1 Constraints on the radial dependence of light deflection 



Mortlock & Turner (2001b) considered a "zero-th" order approach, by recalling 
that the deflection angle in GR is twice the Newtonian value. They simply 
assumed that a similar relation holds for MOND, and derive the "deflection 
law". We note, however, that this approach is rather simplistic, and ad hoc. 
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because the force law used by [Mortlock fc Turner (2001b| ) does not return to 
the Newtonian form on large scales (although with a modified gravitational 
constant), which is required from general considerations (e.g., [Sanders, 1986| ; 
[Walker, 19941 ; [White fc Kochanek, 200l|) . 



The aim of this section is to give an indication of the potential of current data 
sets for the purpose of testing alternative theories of gravity. Along this line 
of reasoning, we first concentrate of the general class of models considered by 
Mortlock fc Turner (2001a| ). One simply starts from the deflection law for a 



point mass. In General Relativity, the deflection law is 

«W = -f , (9) 



where is the Einstein radius of the lens. The Einstein radius depends on 
the mass M of the lens and the angular diameter distances from observer to 
lens (A), observer to source (-D^), and lens to source (As): 



UGM D 



V DiDi, 

Hence, another complication of lensing in alternative theories of gravity arises 
because it is not clear how the angular diameter distances vary with redshift. 
However, this uncertainty should only affect the amplitude of the lensing sig- 
nal, and not its angular dependence. 



Mortlock & Turner (2001a) considered a generic deflection law, which is pa- 
rameterized as 



For 9 <^ 9o, this deflection law reduces to the usual deflection, as is the case 
for ^ = 1. For large 9 the deflection angle decreases oc 9~^. The case of ^ = 0, 
can be considered as the approximate MOND deflection angle. 

It is straightforward to relate the deflection angle to the observable shear. For 
this particular choice of deflection angle we obtain ( [Mortlock fc Turner, 2001"a|) 

Mortlock & Turner (200Ta) compared this model to the observed ensemble 
averaged tangential shear around galaxies in the SDSS (p^'ischer et al., 20lin[). 
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Fig. 7. Left panel: Likelihood contours for ^ = 0, which can be considered as a 
good approximation to the MOND case. The contours indicate the 68.3%, 95.4% 
and 99.7% confidence Umits on two parameters jointly. The allowable range for the 
value of Oq excludes small values: the MOND regime sets in at too large physical 
scales. The best fit value for the Einstein radius corresponds to a mass-to-light ratio 
of M/Lb = I.SMq/Lb©. Right panel: The likelihood contours for the parameters ^ 
and 9q as determined from a maximum likelihood analysis of 16 deg^ of RCS Rq 
band imaging data. The results have been marginalized over the Einstein radius 
^E- As expected, ^ = 1 (GR without dark matter) is ruled out. Also the best fit 
of = 0.1 from Mortlock &: Turner (2001a ) is ruled out with high confidence. In 



addition, the maximum likelihood solution is much worse than the GR + truncated 
halo model presented in section 4, and the current data are sufficient to rule out 
models of gravity with deflection laws given by equation 11 at the lOa confidence 
level. 

They found a best fit value of ^ = —0.1 and = 3.7 arcsecond. As a result, 
[Mortlock fc Turner (200l"a|) concluded the MOND "predictions" are consistent 
with the SDSS measurements. 



However, as we will show now, a maximum likelihood analysis can provide 
much better constraint. As before, we use the measurements from the RCS. 
We now treat the lenses as point masses with the deflection law given by 
equation 11. 

Figure 7 shows the results of a maximum likelihood analysis of 16 deg^ of Rc 
band imaging data from the RCS. The left panel shows the constraints on 
Oq and for ^ = (i.e., the "MOND" approximation). The allowable range 
for the value of 9q excludes small values: the MOND regime sets in at too 
large physical scales. The best fit value for the Einstein radius corresponds to 
a mass-to-light ratio of M/Lb = I.SMq/L^q, which is similar to what one 
would expect from the stars themselves. The right panel of Figure 7 shows 
the constraints on and ^ when we marginalize over 6'e. Both the ^ = 1 



19 



(GR without dark matter) and the ^ = case ("MOND") are excluded by the 
current data. 



Comparison of the maximum hkehhood solution to the best fit "GR + trun- 
cated dark matter halo" model shows that the latter is a much better fit. The 
results from the "truncated dark matter halo" model indicate that the data 
favor a model in which the projected mass distribution has a profile that is 
cx: 1/r on small scales, and steepens on larger scales. This explains why the 
,^ = model results in large values of 6q, and why the data favor a value of 
^ ~ 0.35: the best fit model is a compromise between the 1/r and profiles, 
reflecting the change in slope in the density profile. Comparison of the likeli- 
hoods of the deflection law models given by equation 11 and the "truncated 
halo" model shows that the former models are excluded at the lOo" confidence 
level. 



The problems with the deflection laws considered by [IVLortlock fc 'I'urner 



(2001a|) arise because of the form of the deflection law on large scales. As 
mentioned above, the law of gravity has to return to the 1/6 form on large 
scales, and consequently, such theories should not only specify a scale where 
the effect of alternative gravity becomes important, but also large scale ^out, 
where the theory returns to "normal" gravity. 



To allow for such a scale, we modify the [Mortlock fc Turner (200l"aD deflection 
law (for ^ = 0, i.e., the "MOND" case) to 

66o \ d + dont 



«W = -7if ^ + ^0-7^^ , (13) 



which gives a shear 



m 



9^ 



(14) 



For 9 ^ ^^out this modified deflection law leads to a deflection angle (or shear) 
that is a factor 9out/9o times the GR value. Hence, the ratio 9ont/9o is a measure 
of the "mass descrepancy" . Comparison with the truncated halo model shows 
that both profiles are nearly identical on scales larger than a few 9q, provided 
the truncation parameter s = 6'out- We find that such a model provides a good 
fit to the data. However, we note that the introduction of ^out is rather ad 
hoc. 

The study of the radial dependence of the deflection law is potentially useful 
when it is predicted completely by the theory. The ad hoc approach used here 
can fit the data, and hence gives inconclusive results. As mentioned above, 
the azimuthal variation of the lensing signal around galaxies provides a more 
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robust method. The detection of an anisotropic lensing signal poses a serious 
problem for alternative theories of gravity. 



6 Conclusions 



Weak gravitational lensing has shown tremendous progress in the last few 
years, thanks mainly to the fact that large areas of the sky can now be observed 
with wide field cameras, such as the CFH12k camera on CFHT. In this review 
we have highlighted some of the most recent results. 



The most recent cosmic shear studies ( Paeon et al., 2002| ; [Hoekstra et al. 



2002b ; Refregier et al., 2002 ; Van Waerbeke et al, 200^ ) show excellent agree- 



ment. This is remarkable, given the fact that the data have been taken with 
different instruments, filters, and depths. These most recent results demon- 
strate that weak lensing by large scale structure can play an important role 
in the era of precision cosmology. 

Weak lensing provides a direct measure of the (dark) matter distribution on all 
scales, and hence it is one of the most powerful probes of the relation between 
the galaxies and the underlying mass distribution. As a result, it provides 
a direct way to determine the bias parameter b and the galaxy-mass cross- 
correlation coefficient r as a function of scale, on scales not accessible by other 
techiques. The results obtained by [Hoekstra et al. (2002c| ) suggest significant 
non-linear or stochastic biasing on scales 0.5 — l/i^^ Mpc. These measurements 
provide important observational constraints on models of galaxy formation. 

The lack of visible tracers hampers our knowledge of the gravitational po- 
tential around galaxies at large projected radii. Weak lensing has proven to 
be a powerful tool to probe the extent and shapes of dark matter halos. The 
maximum likelihood analysis of RCS data indicates a clear steepening of the 
mass distribution on scales larger than ~ 250/i~^ kpc ([Hoekstra et al., 2002dD. 



Furthermore, the data allowed for the first detection of the fiattening of dark 
matter halos, with sperical halos excluded at the 99% confidence level. 

Finally weak lensing can be used to constrain alternative models of gravity 
(without dark matter). Ad hoc approaches to study the radial dependence of 
the defiection law remain inconclusive, mainly because of the lack of predic- 
tions. These theories, however, predict an essentially isotropic lensing signal 
around galaxies, irrespective of the actual defiection law. The detection of an 
anisotropic lensing signal around galaxies therefore poses a serious problem 
for theories of alternative gravity, such as MOND. 

The main results discussed in this review were obtained in the last few years. 
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The progress made in such a small amount of time suggests that exciting times 
are ahead, when even larger data sets become available. 



HH thanks the organizers for the invitation to an interesting meeting. 
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